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I.  IMTRODUCTKai 


During  the  past  several  years  Midwest  Research  Institute  has  been 
engaged  in  a  study  of  certain  stress -dependent  aspects  of  iiltrasonic  propaga¬ 
tion  in  solid  materials.  A  thorou^  understaunding  of  these  stress -dependent 
properties  should  greatly  facilitate  the  ultimate  use  of  ultrasonics  in  all 
types  of  stress  measurement  problems  including  the  nondestiructive  measurement 
of  residual  stress. 

The  propagation  velocity  of  ultrasonic  waves  in  solids  varies  with 
stress  in  a  manner  that  depends  on  the  third-order  elastic  constants.  In 
addition,  shear  wave  velocity  depends  upon  the  angle  between  the  polarization 
axis  and  the  stress  direction.  These  properties  can  combine  to  produce  an 
effective  rotation  of  the  polarization  axis  as  a  shear  wave  travels  through 
a  stressed  medium.  Techniques  for  measuring  these  stress-dependent  prop¬ 
erties  have  been  described  in  previous  reportsl>^/*  along  with  e:q>erimental 
results  for  several  different  materials. 

The  stress-dependent  rotation  of  shear  wave  polarization  is  very 
similar  to  certain  photoelastic  effects.  It  also  suffers  many  of  the  same 
limitations.  One  such  limitation  is  that  the  total  rotation  observed  is  the 
integrated  effect  of  travel  through  a  given  thickness.  Thus,  stress  varia¬ 
tions  through  a  given  thickness  may  produce  a  csinceling  effect  and  give 
results  that  are  identical  to  those  obtained  in  a  stress-free  medium.  This 
problem  has  been  solved  in  the  case  of  photoelasticity  through  the  develop¬ 
ment  of  a  three-dimensional  technique?/  which  depends  on  certain  peculiarities 
of  scattering  from  a  polarized  light  beam.  A  somewhat  analogous  technique  is 
now  theoretically  possible  using  ultrasonic  waves  and  the  results  described 
in  this  report. 


II.  STRESS-IHDUCED  BIKEFRINGENCE 


In  the  following  sections  of  this  report  several  areas  of  work  are 
described  relative  to  the  use  of  ultrasonics  in  a  technique  of  three-dimen¬ 
sional  stress  analysis.  As  a  prelude  to  this  study  it  was  felt  that  addi¬ 
tional  information  was  needed  about  the  behavior  of  shear  waves  in  biaxially 
stressed  samples.  In  earlier  work  we  had  limited  our  birefringence  experi¬ 
ments  to  uniaxial  stress  conditions.  Results  of  the  earlier  workiz§/ 
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indicated  that  AV/V  was  linearly  proportional  to  the  uniaxial  stress.  An 
ejqperiment  was  therefore  set  up  to  determine  the  relationship  between  AV/V 
and  P-Q  ,  the  difference  between  the  principal  ntresses  in  a  bleucially 
stressed  specimen. 

The  specimen  shaije  and  loading  conditions  are  shown  in  Figure  1. 
llie  8-in.  dieuneter  disk  was  S/S  in.  thick  and  made  frcrni  G061-T6  aluminum 
plate.  Values  of  AV/V  were  determined  at  several  different  locations  on 
the  face  of  the  plate  using  the  pulse -echo  technique  described  previously. i/ 
The  compressive  load  was  applied  along  a  diameter  of  the  disk  either  parallel 
or  perpendicular  to  the  rolling  direction  of  the  alucdnum  stock.  The  aniso¬ 
tropic  rolling  texture  of  aluminum  produces  a  birefringence  effect  with 
principal  axes  parallel  and  perpendicular  to  the  rolling  direction.  With  the 
loading  parallel  to  one  of  these  principal  axes,  superposition  of  the  original 
state  (anisotropy  due  to  preferred  orientation)  and  the  stress  state  produces 
a  rotation  of  the  resultant  principal  axes  at  all  points  on  the  face  of  the 
disk  except  along  the  two  diameters  that  are  parallel  and  perpendic\ilar  to 
loading.  It  was  along  these  two  (primary)  diameters  that  measurements  of 
AV/V  were  first  made. 

The  value  of  (P-Q)  at  any  point  (x,y)  on  the  disk  can  be  cal¬ 
culated  using  the  expression^/ 


4rF(x^4-y^-r^) 
TT^(x^+y^-r^)  +  4r^x^J 


where  r  is  the  disk  radius  and  F  is  the  diametral  load.  The  results  of 
measurements  at  the  center  of  the  disk  and  one  additional  point  on  each  of 
the  two  primary  diameters  conclusively  showed  that  AV/V  does  vary  linearly 
with  (P-Q)  . 

In  addition  to  measurements  along  the  primary  diameters,  we  made 
several  attempts  to  verify  the  stress  state  at  other  points  on  the  face  of  the 
disk.  As  mentioned  before,  the  principal  stress  axes  at  these  points  rotate 

when  the  load  is  varied.  Thus,  the  value  of  (P-Q)  as  well  as  the  direc¬ 
tions  of  the  princli»l  suces  varies  from  point  to  point  and  from  one  load  to 
another.  The  fact  that  the  axes  at  a  given  location  rotate  as  the  load  is 
varied  makes  it  necessary  to  rotate  the  transducer  in  order  to  evaluate 
AV/V  .  Another  problem  arises  because  the  transducers  cover  an  area  over 
which  considerable  variation  may  occur  in  both  (P-Q)  and  the  principal 
axes'  directions.  These  factors  made  it  difficult  to  obtain  clean  pulse- 
echo  patterns  of  the  type  used  in  accurate  evaluations  of  AV/v  . 
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III.  ORAIK  BOUNDARY  SCAITERINa  FROM  A  SINOIE 
UUmSOmC  EBAM 


The  "scattering"  technique^  of  three-dimensional  photoelasticity 
has  recently  been  refined  to  the  point  of  being  a  very  useful  stress  analysis 
tool.  We  will  not  go  into  the  details  of  this  technique,  but  a  cornerstone 
of  the  whole  procedure  is  the  fact  that  polarized  light  waves  produce  imucLmum 
scattered  radiation  In  directions  perpendicular  to  the  E-fleld  vector,  and 
minimum  scattered  intensity  parallel  to  the  E-field  vector.  It  was  felt  that 
some  of  the  techniques  used  in  three-dimensional  photoelasticity  might  be 
extended  to  the  use  of  polarized  shear  waves  and  thus  permit  stress  analysis 
in  structiiral  metals  as  well  as  transpeurent  models.  It  was  first  necessary 
to  study  beam  scattering  in  polycrystalline  metals  and  determine  if  a  spatial 
distribution  exists  similar  to  that  exhibited  by  polarized  electrooagnetlc 
radiation. 


Our  investigation  of  beam  scattering  has  been  limited  to  geometric 
arrangements  similar  to  that  shown  in  Figure  2.  Initial  e^qteriments  involved 
two  5-mc.  AC-cut  quartz  cr^'stals  mounted  on  a  block  of  aluminum  as  shown  in 
Figure  2(a).  The  arrows  on  the  transducers  denote  polarization  direction. 

The  small  electrical  signals  produced  by  the  receiver  crystal  were  cu^llfled 
and  displayed  on  a  CRT  \diose  sweep  was  Initiated  when  the  transmitter  was 
pulsed.  With  this  arrangement  we  were  able  to  determine  the  time  delay  be¬ 
tween  the  transmitted  pulse  auid  receipt  of  the  scattered  signals.  It  was  not 
surprising  to  find  that  a  period  of  "quiet"  followed  the  sweep  initiation, 
but  it  was  rather  suzprlsing  to  find  that  the  first  significant  signals 
occurred  at  a  time  that  corresponded  to  shear  wave  propagation  along  the  path 
(di+d2)  as  shown  in  Figure  ?(b).  This  fact  was  verified  at  several  values 
of  dj^  and  .  A  few  very  weak  signals  are  sometimes  present  between  tQ 
and  but  they  are  usually  much  weaker  than  those  occurring  at 

t(d^+d2)  '  believe  the  earlier  signals  come  from  beam  spreading,  mode 

conversion,  and  multiple  scatteilng  which  may  permit  a  more  direct  path  of 
travel  between  the  two  transducers. 


The  magnitude  of  the  first  strong  signal  is  very  sensitive  to 
slight  translations  of  either  the  transmitter  or  receiver  crystals.  This  is 
probably  due  to  variations  in  the  grain  structure  that  produces  the  scatter¬ 
ing.  The  lack  of  reproducibility  in  the  received  signal  aiqplltude  made  it 
liiq>oB8lble  to  observe  any  systematic  variation  in  scattered  intensity  as  the 
Ix>larlzation  direction  of  either  crystal  was  varied.  Thus,  we  have  been 
unable  to  definitely  pinpoint  a  spatial  variation  in  the  energy  scattered 
from  a  polarized  shear  wave. 
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Figure  2  -  Schematic  Diagrams  of  (a)  Test  Block  With  Transducers 
and  (h)  Section  Sbodng  Apparent  Path  of  Primary 
Wave  and  Scattered  Signal 


The  study  of  beam  scattering  was  next  extended  to  Include  various 
ccnbinatlons  of  X-cut  and  AC-cut  transmitter  and  receiver  crystals.  In  each 
case,  the  resonant  frequency  of  both  crystals  was  the  same.  The  magnitude  of 
the  received  signal  was  much  greater  when  AC-cut  crystals  were  used  for  both 
transmitter  and  receiver  than  €uiy  other  combination  of  AC-cut  and  X-cut 
crystals . 


In  the  previous  discussions  we  have  been  Interested  only  In  the  In¬ 
tensity  of  the  first  signals  to  reach  the  receiver  crystal.  Subseqiient  sig¬ 
nals  are  numerous  euid  ultimately  non  together  to  form  an  almost  continuous 
"noise"  along  the  CRT  sweep.  Multiple  reflections,  beam  spreading,  and  mode 
conversions  give  rise  to  this  relatively  high  intensity  noise  which  finally 
dies  out  as  the  energy  of  the  primary  wave  is  dissipated. 


IV.  INTERACTION  OF  TWO  UI/TRASQNIC  WAVES 


A.  General  Theory 

The  details  of  the  theoretical  analysis  of  elastic  wave  interaction 
in  nonlinear  solids  were  given  in  a  previous  report^/  suid  therefore  irtll  not 
be  repeated  here.  Let  it  sxiffice  to  say  that  an  interaction  between  two 
waves  is  predicted  and  a  "scattered"  wave  of  combination  (siam  or  difference) 
frequency  theoretically  originates  at  the  region  of  intersection.  The  in¬ 
tensity  of  the  scattered  wave  is  expected  to  be  a  maximum  when  so-called 
resonant  conditions  exist.  Table  1  summarizes  the  -/arious  types  of  inter¬ 
action  predicted  by  theory  emd  gives  esqpresslons  for  (cos  9)  >daere  9  is 
the  angle  between  the  primary  beams  at  resonance. 

In  reference  2,  it  was  pointed  out  that  from  a  quantum  mechanical 
standpoint  the  phonon-phonon  interaction  should  give  rise  only  to  phonons 
with  the  sum  of  the  frequencies  of  the  primary  phonons  due  to  conservation  of 
wave  vector  (ks  »  Tsi-*^)  and  energy  (0)3  »  uu^+tug)  .  It  will  be  shown  in 

the  following  discussion  that  the  classical  treatment  vdiich  predicts  a  scat¬ 
tered  wave  with  the  difference  frequency  (u)2^-a)g)  also  has  its  quantum 

mechanical  analogy.  The  argument  is  based  on  a  paper  by  Slonimskii§/  on 
sound  absorption  in  solids. 

It  seems  obvious  that.  If  it  is  possible  for  two  phonons  to  combine 
to  give  a  third  one  of  frequency  (03  ■  ,  the  reverse  process  can  occur 

as  well,  i.e.,  one  phonon  of  frequency  wj  may  split  into  two  under  the 
condition:  -  »],  +  ((Oj-O)^)  .  If  Is  set  equal  to  mg  ,  the 
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is  the  vel(x:ity  ratio  Ct'^l  • 

!  frequency  of  the  longitudinal  primary  wave  is 


egtilvalence  of  the  two  reXetlone  is  clear.  If  the  vsvelengths  of  the  phonons 
involved  are  sufficiently  greater  than  the  lattice  spacing  so  that  dispersion 
and  Unklapp  processes  need  not  be  considered,  it  is  easily  seen  from  the  con¬ 
servation  lavs  that  only  longitvidlnal  phonons  nay  split  into  two  phonons. 

The  resultant  phonons  nay  both  be  transverse  or  one  may  be  longitudinal  and 
the  other  transverse. 

The  elastic  vaves  are  represented  as  a  sum  of  harmonic  vaves 


T>(#)  -  Z 

a 

with  %(]4,r)  -  ^a(a^e^^  ^+aje  (l) 


vhere  e^  is  a  unit  vector  in  the  direction  of  polarisation,  a^^  is  the 
eunplitude,  and  the  propa{^tlon  vector.  If  represents  the  number  of 
phonons  in  a  particular  mode  a  ,  the  natirlx  elements  of  the  displacement 
operator,  according  to  the  theory  of  the  harmonic  oscillator,  are: 


<Na|a*|N,-l> 


(2) 


■vdiere  m  is  the  mass  of  the  volume  of  interest.  We  are  interested  in  the 
transition  probability  from  a  state  with  nimber  of  phonons  ,  N2  ,  N3  to 
a  state  with  number  of  phonons  Nj^-l  ,  N2+I  ,  Nj+l  .  The  part  of  the  elastic 
energy  expansion  containing  cubic  strain  terms  is  taken  as  the  perturbing 
Hamiltonian  ^  .  It  is  shown!/  that  ^  can  be  evaluated  to  give: 

-  i[2B(a+b)+2Ac]tt^^^pl,f]tS^®^[l^,^]i!^^^[k5,i5  (3) 


for  the  decay  of  a  longitudinal  phonon  into  one  longitudinal  and  one  trans¬ 
verse  phonon,  and 


8 


-  i[2in>f2Ac+(K-2/3>»)d-2^f]u^^^[ki,r]u^^^p2,r]u^®^^^ 


(4) 


for  tine  decay  of  a  longitudinal  jshonon  into  two  transverse  phonons.  A  and 
B  are  third-order  elastic  constants  euid  a  ,  b  ,  c  ,  d  ,  f  are  expressions 
containing  the  wave  vectors  and  polarizations,  i.e.,  a  ■ 

The  matrix  elements  are  then 


<Ni,N2,M3|hJ^^ 


Ni-l,N2+l,N3+i> 


e 


(5) 


for  the  catis  (3)  where  K  =  ki-k2-k3  ,  and  0  ■  a)i-a)2-<03  .  Integrating  over 
the  volume  V  and  calculating  the  transition  probability  per  unit  time 
according  to  perturbation  theory  one  obtains; 


„(1) 


2^  Nj^(N2^1)(H3^1) 
(2m)2  ‘"l“>2'^ 


[2B(a+o)+2Ac]^6(fi)[  d\] 


(6) 


The  volume  integral  goes  rapidly  to  zero  for  K  0  if  V  is  large  compared 
to  the  wavelengths  involved.  When  K  ■  0  ,  we  see  that  is  propor¬ 

tional  to  V®  .  Using  expression  (4)  a  similar  expression  is  obtained  for 
the  transition  probability  of  a  longitudinal  phonon  into  two  transverse 
phonons . 


It  is  now  apparent  how  the  generation  of  an  elastic  wave  of  fre¬ 
quency  (tt)2^-«D2)  >  predicted  classically,  can  also  be  predicted  from  a 
qmntum  mechanical  standpoint.  Furthermore,  the  qusuitum  mechanical  treatment 
allows  an  easy  calculation  of  the  intensities.  Equation  (6)  shows  that  the 
transition  probability  is  proportional  to  the  number  of  ph<»jon8  involved  in 
the  process.  Now,  consider  the  state  |Ni,N2,N^  ,  from  which  the  transi¬ 
tion  to  the  state  |Ni-1,N2+1,N3+1^  will  be  meide.  Note  that  the  transition 
probability  is  enhanced  by  a  large  N2  Just  as  it  is  by  a  large  Ni  .  In 
other  words,  if  two  elastic  waves  with  frequencies  (co^  >  (#2)  intersect 
under  the  right  conditions  as  outlined  above,  the  Interaction  results  in  a 
reduction  of  the  intensity  of  the  wave  .  Hie  amount  is  partly  added  to 


9 


the  Intensity  of  the  wave  <02  end  pertly  gives  rise  to  e  third  veve  of  fre¬ 
quency  .  One  miidit  say  that  the  second  wave  stioiuletes  the  split¬ 

ting  of  the  ]^onons  of  the  first  wave.  This  "maser ''-like  intersustion  may 
have  some  further  significance. 

The  correspondence  between  the  classical  and  qusuitum  mechanical 
treatment  of  two  interacting  elastic  waves  in  a  solid  now  appears  more 
reasonable.  Both  treatments  lead  to  the  same  conditions  for  the  angle  of 
intersection  of  the  two  waves  to  produce  a  third  one  of  either  the  sum  or 
difference  frequency^  and  also  give  the  right  polarizations  of  the  three 
waves  relative  to  each  other.  It  is  expected  further  that,  for  the  same 
intensities  of  the  two  primary  waves,  the  intensity  of  a  scattered  wave  with 
frequency  (a)i-(U2)  than  the  intensity  of  a  scattered  wave  with  fre¬ 
quency  ,  since  in  the  latter  case  two  phonons  combine  to  give  a 

third  one  corresponding  to  the  tnuisltlon  from  a  state  ]%,N2,N^  to  a 
state  |Ni-1,N2-1,N3+1>  . 

B.  Intensity  of  Interaction 

Using  the  techniques  described  in  a  previous  report,  expressions 
have  been  derived  for  the  amplitude  of  scattered  waves  from  all  of  the  pos¬ 
sible  interaction  cases.  These  e;q>ressions  are  given  in  Appendix  A.  Ihe 
most  favorable  conditions  for  detecting  beam  Interaction  should  occur  >dien 
the  anplltude  of  the  scattered  wave  is  maximum.  Obviously  the  maximum  should 
occur  under  conditions  of  resonance,  but  it  should  be  remembered  that  reso¬ 
nance  can  occur  over  a  finite  range  of  (a)  values.  !nie  asplltude  e:q)resslon6 
sho'.fn  in  Appendix  A  were  therefore  programmed  for  an  IBM  1620  computer  and 
numerical  results  were  obtained  for  various  values  of  the  frequency  ratio, 

(a)  ,  and  for  most  of  the  materials  for  >rtiich  third-order  elastic  constants 
eu:%  avallable.^LZ/ 

Figures  3-9  are  plots  of  amplitude  factors  for  four  different 
materials  in  each  of  the  seven  cases  that  yields  a  nonzero  airplltude  expres¬ 
sion.  ISie  antplitude  is  directly  proportional  to  the  factor  y  through  the 
eiqpresslon 


anplitude  s 


(7) 


idiere  and  Xg  are  the  amplitudes  of  the  primary  waves,  fj^  is  the  fre¬ 
quency  in  mc/sec,  V  is  the  vol\mie  of  intersection,  and  R  is  the  dlsteince 
between  the  "point"  of  intersection  suid  the  i>oint  of  observation. 
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0  0.2  0.4  0.6  0.8  1.0 

(a) 

Figure  3  -  AngpUtude  Factors  for  the  Interaction  of  Two  Transverse  Waves 
to  Produce  a  Scattered  Longitudinal  Wave  (Case  I-a) .  Both 
Primary  Waves  Polarized  Peirendicular  to  Plane. 
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(a) 

-  Aaaiplitude  Factors  for  the  Interaction  of  Two  Transverse  Waves 
to  Produce  a  Scattered  Longitudinal  V/ave  (Case  I-B) .  Both 
Primary  Waves  Polarized  Parallel  to  kiic2  Plane. 
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Figure  7  -  Amplitude  Factors  for  the  Interaction  of  a  Longitudinal 
and  a  Transverse  Wave  to  Produce  a  Difference  Fiequency 
Longitudinal  Wave  (Case  IV-A).  The  Prlniary  Transverse 
Wave  is  Polarized  Parallel  to  the 
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(a) 

-  Anqplitude  Factors  for  the  Interaction  of  a  Longitudinal 
and  a  Transverse  Wave  to  Produce  a  Difference  Frequency 
Transverse  Wave  (Case  V-B)  .  The  Primary  Transverse 
Wave  is  Polarized  Parallel  to  the  k2^k2  Plane. 
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Y  has  the  dimensions  of  sec^cm^  in  the  c.g.s.  system  used  here.  We  see 
that  Y  is  free  of  all  of  the  experimental  variables  except  the  frequency 
ratio  (a)  and  the  physical  constamts  of  the  material  under  investigation.  It 
has  been  mentioned  previously?/  that  the  amplitude  expressions  were  derived 
for  the  following  conditions: 


R  »  »  Xg 


where  Xg  is  the  wavelength  of  the  scattered  wave.  These  conditions,  of 
course,  prohibit  R  from  going  to  zero  in  (7)  above. 

In  Figures  3  -  9,  we  have  plotted  the  absolute  values  of  y  .  The 
real  value  of  y  may  be  negative  as  well  as  positive,  but  the  intensity  of 
the  scattered  wave  is  proportional  to  y^  and  therefore  always  positive.  In 
plotting  the  data  we  have  used  dashed  lines  near  those  areas  \rtiere  the  real 
value  of  Y  changes  sign.  The  absolute  value  of  y  always  shows  a  minimum 
near  such  areas  as  it  actmlly  goes  to  zero.  The  exact  locations  of  these 
zero  intensity  points  have  not  been  accurately  determined  because  we  are  more 
interested  in  the  intensity  maxima. 

The  information  contained  in  the  (y  versus  a)  plots  can  be  used  in 
optimizing  the  experimental  variables.  IJje  curves  are  particularly  helpful 
in  locating  relatively  sharp  maxima  such  as  those  shown  in  Figure  6.  In 
general,  we  see  that  the  intensity  of  the  interaction  should  be  larger  in 
polystyrene  than  any  of  the  other  materials  examined  here.  One  might  have 
guessed  this  on  intuitive  grounds  and  we  expect  that  other  "plastics"  would 
probably  comijare  favorably  with  polystyrene.  It  should  be  remembered,  how¬ 
ever,  that  the  calculations  consider  all  materials  as  dissipationless .  Ex¬ 
pression  (7)  indicates  the  scattered  amplitude  has  a  very  strong  dependence 
on  frequency  as  well  as  arqjlltude  of  the  primary  waves.  It  is  possible  (and 
indeed  probable)  that  a  material  with  a  high  y  will  be  particularly  absorp¬ 
tive  and  therefore  attenuate  both  the  primary  and  scattered  waves  much  more 
rapidly  than  materials  with  a  low  y  . 


C.  E3q)erlmental  Details  and  Results 

Figure  10  shows  a  block  diagram  of  a  typical  setup  for  observation 
of  scattered  waves  produced  by  ultrasonic  beam  intewwtion.  The  triangular 
specimen  is  cut  so  that  the  Intersection  angle  9  satisfies  the  resonant 
condition  for  interaction.  Transducers  (A)  and  (B)  are  driven  by  separate  rf 
pulse  generators  and  produce  the  two  primary  waves.  An  Arenberg  pulsed  oscil¬ 
lator  model  PG-650C  is  used  to  drive  one  transducer  and  the  pulsed  oscillator 
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Figure  10  -  Block  Diagram  of  Test  Piece  and  Associated  Electronics 


portion  of  a  Sperry  ultrasonic  attenuation  coogarator  is  used  to  drive  the 
second  transducer.  Transducer  (C)  is  used  as  a  receiver  for  the  scattered 
wave.  At  time  t©  a  signal  from  the  blocking  oscillator  triggers  Pulse 
Generator  No.  1  and  the  oscilloscope  sweep.  A  variable  delay  is  \ised  to 
trigger  Pulse  Generator  No.  2  at  (to+t^)  vhere  t^  is  the  delay.  The 
txinable  amplifier  has  a  bandwidth  of  approximately  1.5  mc/sec  at  the  3  db 
points,  and  operating  into  the  scope  it  produces  a  maximum  deflection  of 
approximately  5  mv/in  at  5  mc/sec.  Quartz  transducers  have  been  used  ex¬ 
clusively;  therefore,  impedance  matching  circuits  have  been  used  between  the 
transducers  and  the  generators  or  ampLtfier. 

Isolation  and  identification  of  the  scattered  signals  are  accom¬ 
plished  in  a  ntmiber  of  ways.  First,  the  time  of  arrival  of  the  scattered 
signal  at  transdicer  (C)  should  correspond  to  acoustic  propeigation  along 
paths  d^  and  d2  .  It  shovild  be  recalled  that  the  velocity  of  propagation 
along  leg  d^^  may  be  different  than  that  along  dg  .  For  example,  the 
triangular  specimen  illustrated  in  Figiire  1  was  made  up  for  the  interaction 
of  two  transverse  waves  of  frequencies  and  (sg  to  produce  a  longi¬ 
tudinal  wave  of  frequency  lOj  »  *  Thus,  the  scattered  signal  should 

arrive  at  (C)  at  time  (tQ+dj^/c^+dg/Cjj)  .  The  time  can  easily  be  determined 
with  the  calibrated  sweep  of  the  oscilloscope.  Transducer  (C)  is,  of  course, 
selected  by  frequency  and  mode  to  match  the  expected  scattered  wave.  Ampli¬ 
fier  tuning  reveals  the  frequency  of  the  detected  signal.  In  addition,  the 
setting  of  t(i  must  be  selected  so  that  both  primary  wave  packets  pass 
through  the  volume  of  interaction  at  approximately  the  same  tine. 

Figure  11  is  an  oscillogram  showing  a  typical  scattered  wave  pulse 
(indicated  by  arrow)  and  Illustrates  the  signal  to  noise  relationship.  The 
signal  amplitude  will  drop  severely  if  (l)  the  amplitude  of  either  primary 
wave  is  reduced,  (2)  the  frequency  of  either  primary  wave  is  changed,  (3)  the 
ampUfier  is  tuned  off  the  expected  0)3  frequency,  or  (4)  tc[  is  varied 
from  the  calculated  value. 

During  the  past  few  months  we  have  experimentally  confirmed  each  of 
the  five  theoretical  Interaction  cases  enumerated  in  Table  1.  The  successful 
detection  of  a  "scattered"  beam  was  first  attained  in  a  triangular  si)ecinen 
of  fused  silica  using  Case  I,  i.e.,  interaction  of  two  transverse  waves  to 
produce  a  sum  frequency  longitudinal  wave.  For  short-hand  purposes  let  us 
refer  to  such  aui  interaction  as  follows:  T((i)2^)+T(»2)  =  L(a)i+u)2)  .  The  actual 
case  studied  in  the  fused  silica  specimen  was  T(5)+T(5)  ■  L(lO).  Theoiy  pre¬ 
dicts  (as  shown  in  Appendix  A)  that  a  scattered  wave  should  be  generated  \dien 
the  two  transverse  waves  are  both  polarized  in  the  (k^,k2)  plane  or  when 
both  are  polarized  perpendicular  to  the  (ki,k2)  Plane*  It  further  predicts 
that  the  signal  amplitude  should  be  zero  when  one  primary  wave  is  polarized 
In  the  (k2^,k2)  plajie  and  the  other  is  polarized  perpendicular  to  the  plane. 
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We  have  experimentally  confirmed  all  these  predictions  \ulng  the  fused  silica 
specimen.  Using  the  same  specimen  ve  have  replaced  the  two  5-mc/sec  crystals 
with  one  3-  and  one  7-mc/sec  crystal.  Althou^  the  sum  fi^quency  is  still 
10  mc/sec  we  have  been  unable  to  detect  a  scattered  signal.  Ilhese  negative 
results  eigree  with  theory,  of  course,  because  the  specimen  was  not  cut  to 
produce  resonant  Interaction  between  beams  of  3  and  7  mc/sec. 

We  have  also  performed  several  experiments  to  determine  the  degree 
of  collimation  exhibited  by  the  scattered  wave.  One  such  experiment  is 
Illustrated  in  Figure  12.  With  the  primary  transducers  arranged  symmetrically 
as  shown  in  Figure  12a,  the  maximvm  scattered  signal  was  detected  with  the 
receiving  transducer  centered  at  position  0  .  When  the  receiving  transducer 
was  moved  1  cm.  either  way  to  position  X  or  Y  ,  the  anplitude  of  the  scat¬ 
tered  signal  fell  to  less  than  one-tenth  the  value  obtained  at  0  .  If  the 
primary  transducers  are  rearranged  as  in  Figure  12b  (and  is  adjusted  to 

the  appropriate  delay)  we  now  find  the  maximum  scattered  signal  at  position 
O'  . 


The  interaction  case,  L(5)+T(5)  =  L(10)  ,  was  studied  in  a  specimen 
of  6061-T6  aluminum  similar  to  that  illustrated  in  Figure  13.  A  10-mc/sec 
longitudinal  wave  was  detected  approximately  18  after  tjj  .  The  cal¬ 

culated  travel  time  along  dj  and  ds  is  shown  to  be  17.8  usee.  Further¬ 
more,  the  scattered  signal  was  observed  only  when  the  delay  time  t^  was 
approximately  3  usee.  This  agrees  well  with  the  calculated  delay  of  3.5 
usee.  In  the  above  case,  it  was  found  that  the  magnitude  of  the  scattered 
signal  was  maximum  vhen  the  polarization  of  the  transverse  wave  was  in  the 
plane  and  zero  when  rotated  to  the  perpendicxilar  position.  This 
result  also  agrees  with  theory  (see  Appendix  A). 

In  the  experiments  Just  described  the  two  primary  waves  were  always 
the  same  frequency.  This  leads  to  the  production  of  a  sum  frequency  equal 
to  2(0(5  ,  a  rather  unfortunate  case  since  the  second  harmonic  of  the  primary 
frequency  is  involved.  In  subs<iquent  experiments  different  primary  frequen¬ 
cies  were  often  used.  Many  of  these  experiments  have  been  performed  on  mag¬ 
nesium  due  to  the  relatively  large  interaction  effects  exhibited  by  that 
material.  The  following  cases  have  all  been  investigated  in  magnesivim. 


T(5)+T(5)  =  L(10) 
L(10)+T(5)  »  T(5) 
L(7)+T(3)  -  L(10) 
L(10)+T(3)  =  L(7) 
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Figure  12  -  Diagram  Shovd.ng  Location  of  Scattered  Signal  for  Different  Transducer  Arrangements 
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Figure  13  -  Interaction  of  Longitudinal  and  Sheau:  Waves  in  Specimen  of  Polycrystalline 


L(15)+L(10)  -  T(5) 
L(15)+T(5)  »  L(10) 
L(10)+T(5)  »  L(15) 
L(5)+T(3)  -  L(8) 
L(5)+L(8)  -  T(3) 
L(8)4^(3)  -  L(5) 


Most  of  the  above  oases  -were  studied  using  disk-like  specimens  similar  to 
that  shovn  in  Figure  14.  Small  flats  are  machined  at  selected  points  on  the 
disk  circumference  to  permit  the  application  of  transducers.  The  flats  are 
perpendicular  to  the  disk  radii;  thus,  interaction  between  the  primary  beams 
always  occurs  at  the  center  of  the  disk.  If  each  flat  has  a  parallel  flat 
on  the  opposite  side  of  the  disk,  as  represented  by  positions  A'  ,  B'  amd 
C  ,  then  several  interaction  cases  can  be  studied.  For  example,  consider 
a  5-mc.  longitudinal  wave,  Ei  ,  intersecting  a  3-mc.  shear  wave,  k2  >  "to 
produce  an  8-mc.  longitudinal  wave,  %  .  The  angle  tpi  must  satisfy  the 
resonant  condition  for  this  interaction  and  the  scattered  wave  travels  in  the 
Is3  direction  as  defined  by  (kj^<-k2)  .  ^  close  e:camination  it  becomes  ob¬ 

vious  that  all  of  the  interaction  cases  listed  in  Table  2  can  be  studied. 


TABIE  2 

INTERACTION  CASES  THAT  CAN  BE  STUDIED 
ON  SPECIMEN  SHOWN  IN  FIGURE  14 


Interaction  Case 

Resonance 

Angle 

Direction  of 
Primary  Waves 

Direction  of 
Scattered 

Wave 

L(5)+T(3)  =  L(8) 

'Pi 

kj^  and  kg 

% 

L(8)+T(3)  a  L(5) 

'P2 

1:3  ^2 

% 

L(8)+L(5)  a  T(3) 

93 

k3  and  k^ 

kg 
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Another  experiment  vas  performed  on  a  maGnesiiun  disk  to  further 
explore  the  collimation  characteristics  of  the  scattered  wave.  An  additional 
flat  was  prepared  on  either  side  of  position  C  such  that  their  normals  were 
tl2*  from  that  represented  by  E3  in  Pigiure  14.  Althou^  a  relatively  strong 
signal  was  observed  at  C  (on  a  6-in.  diameter  disk),  no  signals  at  all  were 
detected  at  X  or  Y.  These  results  definitely  rule  out  the  possibility  that 
the  scattered  signals  we  observe  are  due  to  independent  scattering  of  the  two 
primary  beams  from  grain  boundaries  with  subsequent  "mixing"  of  these  signals 
at  the  receiving  transducer.  Such  scattering  would  not  be  particularly  di¬ 
rective  and  produce  such  a  high  degree  of  collimation.  The  relatively  sharp 
collimation  also  suggests  tliat  the  scattered  wave  is  nearly  plane  and  probably 
suffers  less  than  the  R“^  fall -off  assumed  in  calculating  the  amplitude  fac¬ 
tors  of  Appendix  A. 

The  relative  intensities  of  the  primary  and  scattered  waves  have 
been  examined  but  current  results  are  only  qualitative.  Generally,  the 
voltages  measured  on  the  receiving  transducers  are  40  -  60  db  below  those 
produced  by  priinarj’’  wave  echos  returning  to  the  transmitting  crystals. 

We  have  also  been  interested  in  measuring  changes  in  the  primary 
beam  intensities  caused  by  interaction  effects.  As  discussed  in  Section  IV-A, 
the  intensity  of  both  priinary  beams  might  be  e.'qxjcted  to  decrease  in  those 
interaction  cases  where  the  scattered  wave  is  a  sum  frequency.  In  the  case 
of  difference  frequencies,  hov/ever,  the  intensity  of  the  "splitting"  \nxve 
would  be  expected  to  decrease  while  the  other  would  increase.  Experimental 
efforts  to  observe  these  relatively  snail  changes  in  primary  beam  intensities 
have  thus  far  yielded  negative  results. 

Although  the  intensity  evurves  show  tliat  beam  interaction  should  be 
relativelj'  great  in  polj'sti*rene ,  w  have  been  unable  to  observe  a  scattered 
wave  in  this  material.  These  results  can  probably  be  explained  by  the  large 
ultrasonic  attenuation  e;dilbited  by  polystyrene.  A  high  attenuation  \d.ll 
seriously  reduce  the  intensity  of  the  primary  beams  at  the  point  of  inter¬ 
section  as  well  as  quicl'Jy  attenuate  the  scattered  wave  after  generation. 


V.  VISUALIZATION  OF  ULTRASONIC  WAVE  PACKETS 


A.  Instrumentation 


The  early  attempts  to  detect  scattered  waves  created  by  beam  inter¬ 
actions  were  confined  to  the  use  of  transducers.  These  transducers  are  quite 
sensitive  to  weak  signals,  but  at  high  frequencies  they  are  highly  directive. 
It  was  felt  that  the  negative  results  which  came  first  were  possibDy  due  to 
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slight  errors  in  the  calculated  propagation  direction  of  the  scattered  wave. 
The  transducer  directivity  night  therefore  prevent  detection  of  the  scattered 
wave.  Hie  sesnrch  for  a  less  directive  detector  finally  led  us  into  various 
optical  techniques.  Although  these  techniques  are  limited  to  trsuisparent 
solids,  they  do  have  certain  advantsiges  for  studying  beam  interaction. 

Several  waves  traveling  in  different  directions  csm  be  visualised  at  the 
ssune  time  and  yet  the  sensitivity  of  detecting  one  wa'.'e  in  relation  to 
another  nay  be  adjusted. 

A  simple  schematic  of  oiir  optical  detection  system  Is  shown  In 
Figure  15.  It  Is  basically  a  standard  schlleren  system  except  that  modifi¬ 
cations  have  been  Included  to  permit  the  \tse  of  polarization  techniques. 

When  operated  as  a  simple  schlleren  system,  the  polarolds,  P  and  A  ,  will 
be  removed.  Hie  xenon  flashtubo,  F  ,  produces  a  rather  fuzzy  line  source 
that  is  Imaged  on  an  adjustable  slit,  S  .  The  light  passing  through  the 
slit  is  collected  by  the  lens,  ,  and  converted  into  a  parallel  beam.  TOie 
bean  passes  throiagh  the  test  block,  M  ,  and  then  into  lens,  ,  which 
produces  a  real  image  of  the  slit  source  S  at  S'  .  A  "stop"  is  normally 
positioned  at  S'  to  intercept  the  image  of  the  slit.  However,  ultrasonic 
waves  in  the  test  block  my  produce  localized  gradients  in  the  index  of 
refraction.  These  gradients,  ^en  normal  to  the  optic  axis,  cause  some  of 
the  light  passing  through  M  to  be  scattered  out  of  the  parallel  beam. 

Ihese  rays  pass  by  the  stop  as  illustrated  by  the  dashed  lines  and  are  focused 
onto  a  ground  glass  plate  or  film  at  G  .  The  regions  in  the  test  block  where 
refractive  index  gradients  exist  will  therefore  appear  as  bright  areas  on  a 
dark  field.  When  a  slit  source  is  used,  maximum  sensitivity  of  detecting 
waves  is  achieved  •vdien  the  wave  front  is  peuallel  to  the  slit.  This  factor 
can  be  used  in  interaction  experiments  to  accentuate  the  inage  of  the  scat¬ 
tered  wave  over  that  of  the  primary  waves. 

When  the  polarizer,  P  ,  smd  analyzer,  A  ,  are  used  as  shown  in 
Figure  15,  the  stop  S'  is  not  necesoary.  If  P  and  A  are  "crossed,"  the 
image  at  G  will  generally  be  a  uniform  dark  field.  An  ultrasonic  wave  in 
the  test  block  M  ,  however,  may  "depolarize"  the  parallel  besun  and  thus 
permit  some  light  to  pass  throxagh  the  analyzer.  IRje  location  of  waves  will 
appear  as  bright  areas  on  a  dark  field.  Wteeks®/  has  shown  that  the  emergent 
light  intensity  is  greatest  for  a  longitudinal  mode  >dien  the  plane  of  polari¬ 
zation  of  the  incident  light  is  inclined  at  45*  to  the  direction  of  sonic 
propagation.  The  modulation  vanishes  when  the  plane  of  polarization  is  paral¬ 
lel  or  perpendicular  to  the  direction  of  sonic  propagation.  Just  the  opposite 
is  true  for  transverse  waves. 

In  an  experiment  in  which  we  wish  to  detect  a  relatively  weak  wave 
generated  at  the  intersection  of  two  stronger  waves,  the  use  of  continuously 
generated  >aves  is  prohibited.  It  would  be  very  difficult  indeed  to  absorb 
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Figure  15  -  Optical  System  for  the  Observation  of  Ultrasonic  Wave  Packets 


all  of  the  energy  in  the  primary  'uavee  before  it  is  repeatedly  reflected  at 
surface  boundaries,  nw  reflected  vaves  vould  ultimately  produce  index  gradi¬ 
ents  over  the  entire  specimen  that  vould  probably  be  larger  than  those  pro¬ 
duced  by  the  scattered  vsve.  This  difficulty  can  be  largely  overcome  by 
using  pulsed  vJLtrasonlc  vaves  and  a  pulsed  light  source  for  the  optical  sys¬ 
tem. 


The  blocking  oscillator  controls  the  repetition  rate  of  both  the 
rf  pulse  generator  and  the  pulsing  of  the  light  source.  However,  a  variable 
delay  between  the  blocking  oscillator  cmd  the  flashtube  power  supply  nakes 
it  possible  to  observe  the  wave  packet  after  it  travels  a  given  distance  into 
the  test  specimen.  With  a  given  delay  between  the  pulsing  of  the  ultrasonic 
generator  and  the  flashtube,  one  can  "stop”  the  traveling  packets  at  a  given 
position  and  visually  observe  or  photograph  the  rcs\altant  patterns.  A 
General  Radio  "Strobotac"  was  used  as  the  flashtube  power  supply.  The  ra-31 
xenon  flashtube  (Edgerton,  Oermeshausen,  and  Grier,  Inc.)  has  an  optical 
window  which  facilitates  imaging  of  the  rather  fuzzy  line  source,  and  it  can 
be  flashed  at  rather  high  repetition  rates. 


B.  Photographic  Results 

Plgvu^  16  is  a  sequence  of  photographs  showing  a  7 -me.  shear  wave 
packet  at  various  stages  of  travel  in  a  glass  specimen.  Standing  waves  are 
clear  in  (c)  due  to  superposition  of  the  Incident  and  reflected  traveling 
waves.  Side  lobes  have  developed  in  (d)  as  the  packet  travels  back  toward 
the  transducer.  Figure  17  shows  two  7 -me.  shear  beams  intersecting  at  ap¬ 
proximately  90*.  Interference  between  the  two  beams  is  easily  seen  even  in 
areas  outside  the  region  of  main  lobe  overlap.  The  duration  of  the  light 
pulse  in  both  figures  was  0.8  asec  corresponding  to  a  shear  wave  travel 
distance  of  approximately  2  mm. 

During  the  past  year,  beam  interactions  were  studied  simultaneously 
using  both  transducers  and  the  optical  system.  Successful  detection  of  the 
sought-after  beams  using  transducer  techniques  lessened  our  immediate  interest 
in  optical  techniques.  Further  experimentation  was  therefore  delayed  \mtil 
more  infomatlon  on  the  interaction  process  could  be  collected.  The  prelimi¬ 
nary  results  obtained  with  transducers  now  suggest  that  under  ideal  conditions 
one  shoijld  be  able  to  detect  the  scattered  waves  with  an  optical  system. 
Further  work  along  this  line  is  warwmted  because  of  the  versatility  of  the 
optical  system  in  studying  beam  spreading  and  structxu^  variations  as  the 
scattered  waves  propagate  away  from  the  point  of  origin. 
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Figure  16  -  Sequence  of  Photographs  Shomng  Ultrasonic  Pulse  Traveling 

Through  Glass  Specimen 
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Figure  17  -  Intersection  of  Two  7-mc.  Shear  Wave  Packets  in  a  Glass  Specimen 
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VI.  THBBE-DlHEIMSIOilAI.  STRESS  ANAUEIS  WITH  UUFRASONICS 


The  Bucceesful  generation  and  detection  of  scattered  beans  of 
reasonable  intensity  provides  us  with  a  theoreticeOIy  possible  technique  of 
performing  three-dimensioned  stress  measurements.  Note  in  Table  1  that  a 
transverse  scattered  wave  can  be  generated  in  either  case  II  or  case  III. 

For  sintplicity,  let  us  consider  only  case  II  \Aere  the  interewjtion  of  two 
longitudinal  waves  may  prod'use  a  transverse  wave  of  difference  frequency.  It 
can  be  shown  theoretically  (and  confirmed  experimentally)  that  the  polariza¬ 
tion  of  the  transverse  wave  lies  in  the  plane  defined  by  the  propagation  vec¬ 
tors  of  the  two  prinary  waves.  We  therefore  have  a  potential  technique  of 
generating  transverse  waves  of.  known  polarization  at  any  subsurface  volume 
element. 


Incorporation  of  beam  Interaction  techniques  into  a  possible  method 
of  three-dimensional  stress  analysis  is  illustrated  in  Figure  16.  Primary 
beams  and  k2  are  first  crossed  at  position  (a)  to  produce  a  transverse 

wave  k3  .  With  a  detector  at  position  (c),  the  total  angular  rotation  o 
of  the  polarization  axis,  can  be  determined  for  pre^gation  between  (a)  and 
(c).  The  primsjry  waves  can  then  be  crossed  at  (b)  and  the  rotation  3  deter¬ 
mined  for  path.  (be).  The  difference  between  the  secondary  principal  stresses 
(in  a  plane  perpendlcxilar  to  its  )  can  then  be  approximated  for  a  point  mid¬ 
way  between  (a)  and  (b).  The  approximation  should  take  the  form 

(p-q)  =  K(a-3)f 

where  p  and  q  are  the  secondary  principal  stresses,  K  is  the  stress- 
rotatlon  constant  for  the  material  in  question  and  I  is  the  distance  between 
(a)  and  (b).  Since  the  direction  %  and  its  polarization  can  be  varied 
through  control  of  Iq.  ^2  >  i'*'  shoiild  be  possible  to  determine  the  di¬ 

rections  of  p  and  q  and  the  value  (p-q)  for  several  different  £3  di¬ 
rections  passing  through  a  common  point. 

There  are  many  practical  problems  to  be  solved  before  the  above  pro- 
cedxrre  will  yield  much  useful  information.  On  certain  geometric  shapes  these 
problems  will  be  much  more  severe  than  others.  The  difficulty  of  varying  the 
direction  of  and  Eg  acciuately  controlling  the  intersection  point 
is  indeed  great,  especially  in  irregularly  shaped  specimens.  However,  liquid 
coupling  techniques  may  prove  satisfactory  in  memy  cases  and  thus  Increase 
versatility  over  contact  bonding.  At  any  rate,  the  potential  applications  of 
beam  scattering  techniques  in  the  measurement  of  stress  and  other  subsurface 
properties  are  numerous  and  deserve  considerable  attention.  Additional  work 
should  be  directed  toward  the  evalmtlon  of  all  influencing  factors  as  well  as 
the  development  of  techniques  that  will  facilitate  the  application  to 
practical  problems. 
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Figure  16  -  Schetfatlc  Diagram  of  Potential  Stress  Analysis  Technique  Using 
t  Beam  Interaction  to  Produce  Sheau:  Waves 
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APPENDIX  A 


AMPUTODE  EXPRESSIONS  FOR  SCATTERED  WAVES 

In  the  ainplltude  expressions,  the  terms  are  defined  as  foUovs: 
X;j^  ,  Xg  -  amplitudes  of  primary  vaves 
V  -  volume  of  intersection 
(I)  -  angular  frequency 
p  -  density 

R  -  distance  betvreen  point  of  Intersecticai  and  point  of  observation 
Ci  -  velocity  of  shear  wave 

Cjj  -  velocity  of  longitudinal  wave 

c+ 

c  -  ratio  of  velocities 

c/ 

^  -  shear  modtilus 
K  -  tjuUt  modulus 

a  -  ratio  of  priznory  frequencies  _ 

<"1 

A  ,  B  -  third  order  elastic  constants 
0  -  angle  between  primary  waves 
%  ’  ^  “  propagation  vectors  for  primary  waves 
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Case  I  -  Two  Transverse  Waves  Interacting  to  Produce  a  Scattered  Longitudinal 
lifcive  of  Prequency  (ars.+tt^) 

(a)  Both  transverse  vaves  polarized  perpendicular  to  the 

plane. 


-  ^K+  5  j  |^c^(3c^-l)(a^+a)+c®(c^“l)(a^+l)J 


16  npR  c; 


(B)  Both  transverse  vaves  polarized  in  the  plEuie. 


antp 


-XiXgVuj 
8  TTpR  Cjj 


ao^(l+a)  Ij-  ^Pi+  ^  +B^  +  cos^ 


U+A+K>2B^ 


(C)  One  transverse  vave  polarized  in 

dicxilar  to 


?3ati=  and  one  perpen- 


aap  a  0 


Case  II  -  Two  Longitudinal  Waves  Interacting  to  Produce  a  Scattered  Transverse 
Wave  of  Frequency  (lug-w^^) 


anp 


-X^XgVi 


8  TTpR  c| 


U+A+K+2^  cos  0  [^-(l+P)cos  0(2a2+2aP)+a2(i+p)2+(a2+p)‘^J 


vhere 


2a^+a^c^2a^c^-a^~6a^2a+c^-l 

2(a-l)^ 
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Case  III  -  One  Transverse  Wave  and  One  Longitudinal  Wave  Interacting  to 
Produce  A  Scattered  Longitudinal  Vteve  of  Frequency  (u)i+a>2) 


(a)  If  the  transverse  wave  is  polarized  in  the  plane. 


airp 


8  npR  4  cl 


a(l-co6^ 


1+a 


{34+A+B)ac+  ^K+  I  +  j  +2B^  c^Cl+a)^  cos  0 


+  V-  ^  A+2BH^  c^  cos  0  +  ^u+  (2a^+4ac  cos  0)  cos  0 


plane. 


(B)  If  the  transverse  wave  is  polarized  perpendicular  to  the  kik2 


emp  «  0 


Case  IV  -  One  Transverse  Wave  and  One  Longltxidinal  Vfeive  Interacting  to  Produce 
A  Scattered  Longitudinal  Wave  of  Frequency  (uJi-U)2) 


(a)  If  the  transverse  wave  is  polarized  in  the  k^kg  plane. 


(B)  If  the  transverse  wave  is  polarized  perpendicular  to  the  k]^^ 

plane. 


amp  s  0 
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Case  V  -  One  TJansverae  Wave  and  One  Longitudinal  Wave  Ii  tertictlng  to  Produce 
A  Scattered  Transverse  Wave  of  Pnsq,uency  (n^-ojg) 


(A)  If  the  prinary  transverse  vave  is  polarized  perpendicular  to 
the  plane . 


amp 


XiXgVd)?  \  1 

16  npR  c^  aj 


U+  ^(4a^+c^-2c^+l+4ac^-4a) 


i  |i+  ^  (2ac®+c^-c®) 


(B)  If  the  primary  transverse  ■wave  is  polarized  in  the  ^2^  plane. 


an^ 


(o^V 


4  npR  c|  Cl 


+  2Mil^  S  COB  0 


vdiere 


M, 


Dj^-(D2"Di)®c  CO® 


M2  = 


(1-a) 


Di  =  - - - j  . 

2  c|(l-cos^  cp) 

"  i 

2c  c^(l-cos2  0) 


2  ^  ja^(2  cos*^  0-1) -ac  cos^  0j  +  p-K-]^ac  cos  0! 


-2(^U+  l^a  cos^  04-  ^Kt-  I  +  I  +2B^  0  cos^  0+(t4-B)c 


Always  use  the  positive  root  of  (l-cos^  0)^ 
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